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Abstract. A two-lane exclusion process is studied where particles move in the two 
lanes in opposite directions and are able to change lanes. The focus is on the 
steady state behavior in situations where a positive current is constrained to an 
extended subsystem (either by appropriate boundary conditions or by the embedding 
environment) where, in the absence of the constraint, the current would be negative. 
We have found two qualitatively different types of steady states and formulated the 
conditions of them in terms of the transition rates. In the first type of steady state, 
a localized cluster of particles forms with an anti-shock located in the subsystem and 
the current vanishes exponentially with the extension of the subsystem. This behavior 
is analogous to that of the one-lane partially asymmetric simple exclusion process, 
and can be realized e.g. when the local drive is induced by making the jump rates 
in two lanes unequal. In the second type of steady state, which is realized e.g. if 
the local drive is induced purely by the bias in the lane change rates, and which has 
thus no counterpart in the one-lane model, a delocalized cluster of particles forms 
which performs a diffusive motion as a whole and, as a consequence, the current 
vanishes inversely proportionally to the extension of the subsystem. The model is 
also studied in the presence of quenched disordered, where, in case of delocalization, 
phenomenological considerations predict anomalously slow, logarithmic decay of the 
current with the system size in contrast with the usual power-law. 
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1. Introduction 

Transport in cells is an intensively studied field of cell biology as it is important 
for cellular organization and function, moreover, malfunction of transport processes 
may lead to severe diseases [H El [3]. Intracellular transport is realized by a class of 
proteins called molecular motors which are able to convert chemical energy released form 
hydrolysis of ATP into mechanical work and thereby perform directed movement on the 
filament network and carry different cargoes from one part of the cell to the other [2]. 
Although basic features of this system are well known there are still important details 
which have not cleared up yet and, for instance, concurring models of bidirectional 
transport exist currently [3]. 

This important issue has motivated much theoretical work on special driven lattice 
gas models which have been used as minimal models of cooperative behavior of molecular 
motors [1]. The cornerstone of these investigations is the asymmetric simple exclusion 
process (ASEP) [5], E] in which particles sitting on a one- dimensional lattice jump 
stochastically to a neighboring site provided the latter is empty. Due to its simple 
formulation, this process has become a paradigmatic model of interacting many-particle 
systems which have nonequilibrium steady-states [3 El E] . The system, when coupled 
to particle reservoirs at the ends, shows boundary induced phase transitions [10] and 
the phase diagram is exactly known [HI [121 Il3] • 

The ASEP has been generalized in many different directions in order to take into 
account more details of real transport systems. Such generalized models are the two 
or multilane models which consist of two or more ASEPs coupled in parallel, where 
particles are able to jump to other lanes. Formally, on the grounds of the directionality 
of the motion in the lanes, these models can be divided in three classes. First, particles 
move totally asymmetrically (unidirectionally) within lanes and the direction of motion 
is the same in all lanes [m [151 IlHl EZl [HI [IS]. This generahzation is inspired by the 
circumstance that the microtubule filament that serves as a track for certain molecular 
motors is not strictly one-dimensional but consists of many proto-filaments. Such models 
have also been investigated in the context of vehicular traffic on multi-lane roads; for 
a review, see e.g. Ref. [20]. Second, two-lane models with unidirectional motion in 
one lane and symmetric (diffusive) motion in the other one have also been introduced 
[211 [22]. These models have been intended to take into consideration that real molecular 
motors can unbind from the filament, diffuse in the surrounding fluid and rebind again to 
the filament. The unbinding and rebinding has also been taken into account by simpler 
models where the ASEP is extended by creation and annihilation of particles [23j, as 
well as by more realistic models where the ASEP is coupled to a finite compartment in 
which the motion is diffusive [2^ [25l[26l[27] . Third, a two- lane model with unidirectional 
motion in the lanes but with opposite orientation of the lanes has also been studied 
[28l [29] . This arrangement can be thought of as a simplified model of motors moving 
along two oppositely oriented filament tracks which are located in a tube-like narrow 
compartment. The motion of a single particle in this environment has been found to 
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have an enhanced diffusion coefficient compared to that of the symmetric random walk 
[29] . This type of active diffusion can be reahzed experimentally and it is also present in 
vivo f3]. An alternative interpretation of this model is that a single filament is placed in 
a narrow compartment where there is a steady convection of the medium in the opposite 
direction compared to the polarity of the filament. 

The theoretical studies have concluded that, in all three cases, the steady-state 
behavior is richer than that of the one-lane ASEP, furthermore, interesting phenomena 
take place in such systems, for instance, synchronization or localization of density shocks. 
The latter can be observed experimentally in the traffic of kinesin motors in accordance 
with theoretical predictions [30, [311 [32] . 

This work wishes to contribute to the studies on the less examined bidirectional 
models of the third class. Earlier studies have concentrated either on the one particle 
dynamics [29] or on the case of weak coupling between the channels, which can be 
treated by means of a hydrodynamic description [28] . In the present work, the case will 
be considered where the inter-lane jump rates are comparable with the intra- lane jump 
rates. A special property of this model is that, under certain conditions, the stationary 
current can be reversed when the global density of particles is varied. In other words, 
if particles are present in the system with a sufficiently high density they are able 
to flow in the opposite direction as a single particle. We shall give the condition of 
this phenomenon in terms of the transition rates of the process. Subsequently, we shall 
analyze the quantitative consequences of this property in arrangements where the system 
contains a finite region where the local bias is opposite to the global one. It will turn out 
that the steady state at such an inclusion with reverse bias, or "barrier" is qualitatively 
different from that of the analogous one-lane partially asymmetric simple exclusion 
process in the reverse bias regime. In the latter model, a shock forms in the middle of 
the barrier and the stationary current tends to zero exponentially with the size of the 
barrier [3ll [331 [I3j . As opposed to this, in the two-lane model to be studied, a shock still 
appears but, under certain conditions, it may be delocalized, which leads to that the 
current vanishes much slower: it is inversely proportional to the size of the barrier. A 
direct consequence of this slow mode is that, in case of quenched disordered lane-change 
rates, the current J vanishes with the size of the system L as J{L) ~ (lnL)~^ in the 
presence of a global bias, which is slower than the usual power-law decay characteristic 
of disordered partially asymmetric one-lane models [351 ESI [371 [381 139] . 

The paper is organized as follows. The model will be defined in Sec. 2. The process 
in the presence of a single particle is studied in Sec. 3. The periodic system and the 
phases determined by the bulk rates are discussed in Sec. 4. In Sec. 5, the phase 
diagram of the open system is determined. Sec. 6 is devoted to heterogeneous models 
including one-barrier systems and random disorder. Some calculations are presented in 
the Appendix and the results are summarized in Sec. 7. 
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The mode is defined on two parallel one dimensional lattices with L sites, denoted 
by A and B, the sites of which are either empty or occupied by a particle. We shall 
study periodic chains (or rings), where site 1 and L are neighbors, as well as open 
chains. The configuration of the system is specified by the set of occupation numbers 
nf'^, I = 1,2, ... ,L, which are zero (one) for empty (occupied) sites. On the space 
of configurations, a continuous-time Markov process is defined by the rates of allowed 
transitions, are the followings for periodic chains. On site / of lane A{B), a particle 
attempts to jump to the adjacent site on its right(left) hand side with rate pi{qi) and 
the attempt is successful if the target site is empty. Furthermore, a particle residing at 
site I of chain A{B) jumps to site I of chain B{A) with rate ui (vi), provided the target 
site is empty. The model is illustrated in Fig. [H 

For open systems, the above process is modified at surface sites 1 and L as follows. 
On the first site of lane A particles are injected with rate a, provided this site is empty. 
On site L of lane A particles are deleted with rate 7, whereas on the first site of lane 
B they are deleted with rate (3. In case of open systems, we restrict ourselves to a 
homogeneous bulk, i.e. pi = p, qi = q, ui = u, vi = v, see Fig. [21 Moreover, we shall 
focus on the reverse bias regime, where a single particle is driven to the left in the bulk 
and are mainly interested in under which conditions a positive current can be forced 
to the system by injecting particles at site 1. For sake of simplicity, particles are not 
injected at site L and the exit rate 7 at that site is set to 1. 
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Figure 1. Schematic view of tlie model in the bulk. Arrows indicate the allowed 
transitions. 



3. One-particle dynamics 

Let us examine the motion of a single particle in the system, a situation which will be 
relevant for the analysis of the stability of shocks. Although, detailed-balance does not 
hold in this system as the motion of the particle is unidirectional within the lanes, we 
shall show that an effective potential can still be defined. Let us consider the steady 
state of a particle on a finite ring. The probability that the particle is found at site i of 
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Figure 2. Schematic view of the open system. Arrows indicate the allowed transitions. 

lane A{B) is denoted by Pi{TTi). The probability current J flowing parallel to the lanes 
(see Fig. [3l) can be written as 

J = Pi-iPi-i - qiiTi. (1) 
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Figure 3. Left: stationary probabilities on a finite ring with a single particle. Right: 
partially asymmetric one-lane model obtained by the elimination of lane B. The 
stationary current across the dashed line is denoted by J. 



Conservation of probability at site i of lane A and B imply 

Pi-ipi-i + Villi - (pi + Ui)pi = (2) 
qi+iTTi+i + UiPi - {qi + Vi)'Ki = 0. (3) 

Let us assume that Vi > for i = 1,2, . . . , L. Expressing vTj from Eq. ([2]) and substituting 
into Eq. ([1]) we obtain 

J = Pi-iPi~i - qipi, (4) 

where = Pi-i{qi + Vi)/vi and qi = qi{pi + Ui)/vi. Formally, this equation is identical 
to that describing the stationary current in a one-lane partially asymmetric model 
with forward and backward rates pi and qi, see Fig. [31 This correspondence with a 
one-lane model is not perfect as the occupations pi alone are not normalized to one, 
J2i=iPi < Y.i=i{Pi + TTj) = 1. Nevertheless, the current differs from that of the 
equivalent one-lane model only by a positive 0(1) factor. This loose formulation of the 
correspondence will be sufficient for latter reasoning as the only information we need is 
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the sign of the current. For the equivalent one-lane model a potential Uf at site i can 
be defined by AUi = Uf — Uf^^ = ln(gj/pj_i). This reads in terms of the rates of the 
two-lane model: 

AU.^Uf-Uf_, = \n ^'(^- + "-) (5) 

Pi-i{qi + Vi) 

In fact, the potential itself exists even in the case when not all fj are positive. This 
effective potential can be obtained in a more direct way by considering the steady state 
of a single particle in a closed system, i.e. open chain with a = /3 = 7 = 0. In this case 
J = and Eqs. ([HIS]) can be equivalently written as: 

qi{pi + Ui) 

Pi-i Pi = ] , r (6) 

Pi-i[qi + Vi) 

^.M = (8) 

qi + Vi 

From these relations, one can see that the steady state is identical to that of an 
equilibrium system, where detailed balance is satisfied and the potential is given by 
Eq. ([5]) in lane A and by Uj^ — = ln[(pj + Ui) / {qi + Vi)] in lane B. We thus conclude 
that the motion of a single particle can be effectively described by a one- dimensional 
motion in the potential Uf". 

Next, we examine the dynamics of a single vacancy in the system when all but 
one of the sites are occupied. It is easy to see that the dynamics of a lonely vacancy 
differs from that of a lonely particle in that the direction of transitions (i.e. the arrows 
in Fig. [3]) are reversed. This is equivalent with a single particle dynamics analyzed 
above, however, in a modified environment where the rates and of the original 
environment are interchanged for all i. The effective potential associated with a 
single vacancy is thus given by 

ATT -TtB jfB , Pi-ljqi+l , . 

qi{Pi-i +Vi) 

In a system which is homogeneous in the bulk a single particle is driven from left to 
right if AU < 0, whereas it is driven from right to left if AU > 0. In the latter case we 
speak of reverse bias and in the followings, we shall focus on this regime. The condition 
Af/ > of reverse bias in the bulk reads in terms of jump rates as : 

q/p > v/u. (10) 

Note that the potential difference experienced by a vacancy can have any sign in the 
reverse bias regime. Namely, it is driven to the right if AU < 0, or equivalently if 
q/p > u/v while it is driven to the left if q/p < u/v. This contrasts with the one-lane 
partially asymmetric exclusion process, where AUi = —AUi always holds. 



Dynamics at barriers in bidirectional two-lane exclusion processes 



7 



4. Periodic system 



Let us consider a homogeneous, periodic system witli 2L sites and particles. It is 
straightforward to check that the probabihty of configurations {nf, nf} is of factorized 
form in the steady state: 

L 



Pi{nt,nr}) 



TT 



where S{i,j) is the Kronecker symbol and the factor Z = I]{„a „s} P({nf , nf})6{Y,i=iinf+ 
nf),N) ensures normalization. The parameters p and n in Eq. ffTTl) lie in the interval 
(0, 1) and satisfy the relation 

up{l-7T) =vn{l- p). (12) 

Keeping the density of particles g = N/{2L) constant and performing the 
thermodynamic limit L —>■ oo, the factorized form in Eq. (ITTl) results in a simple 
form of the current: J^o = pp(l — p) — — ti")- Here, p and vr are fixed by the density 
through g = (p + 7r)/2 and are interpreted as the particle densities in lane A and B, 
respectively. The current, as a function of p, reads as 

qv/u 



Jooip) = P(l -P) 



P- 



(13) 



p+{v/u){l-p)\ 

The fundamental diagram of the model, i.e. the current plotted against the density g, 
can be seen in Fig. HI Depending on the bulk rates, one can distinguish between two 




Figure 4. The current as a function of the particle density g for two different set of 
bulk rates. 



domains of the parameter space in the reverse bias regime, which are characterized by 
quahtatively different shape of the curve J{g). It is straightforward to show that if 

q/p > u/v, (14) 

the fundamental diagram has a single extremum (minimum) and the current is negative 
for any density < ^ < 1. If condition ( |T^ holds, we shall speak of a strong bias (SB). 



Dynamics at barriers in bidirectional two-lane exclusion processes 8 

Note that inequalities fllOl) and fll4p imply that q > p is a necessary condition of strong 
bias. If, however, 

q/p < u/v (15) 

holds, there exists a non-trivial density Qq in the interval (0,1), at which the current 
is zero, i.e. Joo{qo) = 0. The corresponding densities in lane A and B can be easily 
calculated: 




In this case, the fundamental diagram has two extrema (a maximum and a minimum) 
and, for high enough densities {qq < g < I) , the current is positive. If inequality ( ITSi) is 
satisfied, we shall speak of a weak bias (WB). Inequalities (fTOj) and (fT5|) imply that the 
necessary condition of weak bias is v < u. Recalling the dynamics of a single vacancy 
discussed in the previous section, we can see that, in the strong bias regime, the vacancy 
is driven from left to right (AC/ < 0), whereas, in the weak bias regime, it is driven 
from right to left {AU > 0). The connection between the direction of motion of a single 
vacancy and the shape of the fundamental diagram can be immediately seen. If the 
vacancy moves from left to right {AU < 0), the current is negative for finite L and 
N = 2L — 1. Therefore the slope of the curve Joo(f?) must be negative a.t g = 1, i.e. 
lim^^i '^'^'^^^^ < 0. If the vacancy moves from right to left (AU > 0), the current is 
positive and lim^^i '^'^^^^^ > 0. 

In the borderline case q/p = u/v, the motion of a vacancy is diffusive {AU = 0) 
and lim^^i '^'^'^^^^ = 0. Apart from this anomaly at ^ = 1, the shape of the fundamental 
diagram is similar to that in the strong bias regime. 



5. Open system 

In this section, we shall investigate the steady state of an open system in the reverse 
bias regime, i.e. when the bulk rates satisfy inequality (fTOl) . 

First, we recall the steady state of the one-lane partially asymmetric simple 
exclusion process in the reverse bias regime {q > p) [SH [33], |13] . For finite injection rate 
and zero exit rate at site 1, an anti-shock (or domain wall) develops in the density profile, 
where the density changes rapidly with the position (see Fig. [5]). As a consequence of 
particle-hole symmetry of this model, the anti-shock is located in the middle of the 
system. With the distance measured from the middle of the system, the density tends 
exponentially to one(zero) on the left(right) hand side of the anti-shock, and the shock 
region, where the density differs significantly from the bulk value has a characteristic 
width ^ ~ 1/AU = l/ln{q/p). The stationary current decreases with the system size 
as J{L) ~ Q-^^u/"^ and tends to zero in the limit L oo [13]. 

Returning to the two lane model, let us investigate the possibility of an anti- 
shock type steady state. In the low density domain (i.e. on the right hand side of 
the domain wall) and far away from the shock region, particles meet each other very 
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Figure 5. Schematic density profile of the one-lane partially asymmetric simple 
exclusion process in the reverse bias regime q > p. 

rarely. Therefore, they can be regarded as freely moving particles. The same is true 
for vacancies in the high density domain (i.e. on the left hand side of the anti-shock). 
Within the frame of this phenomenological description, we can formulate necessary 
conditions of the stability of the domain wall. First, particles which have come off from 
the high density domain and penetrated into the low density domain must be driven 
back to the high density domain. This condition is equivalent to AU > or, in terms 
of jump rates, to inequality ffTOl) . so it is always satisfied in the reverse bias regime. 
Second, vacancies which penetrate into the high density domain must be driven back to 
the low density domain. This condition is equivalent to AU < or, in terms of jump 
rates, to inequality (HM . We thus conclude that an anti-shock type steady state in the 
reverse bias regime is possible only in the case of a strong bias. This analysis shows 
that the steady state of the model for strong bias and for weak bias are much different; 
therefore the two cases will be discussed separately. 

5.1. Strong bias 

First, we assume that the bias is strong, i.e. the bulk rates satisfy inequality (fHI) . The 
exit rate 7 is set to 1 and we are interested in the steady state for given boundary rates 
a > and (3. Although, we could not find the exact form of the steady state, the shape of 
the density profile on a macroscopic scale and the leading size dependence of the current 
can be obtained by means of phenomenological arguments. As the localized domain wall 
is stable for strong bias we expect the density profile to be qualitatively similar to that 
of the one-lane partially asymmetric simple exclusion process. For (3 = and for any 
a > 0, an anti-shock develops in the interior of the system, which separates a high 
density domain from a low density domain. The center of the shock region is located 
at some < I* < L and the density tends exponentially to zero(one) in the low(high) 
density domain on a characteristic scale ^ ~ 1/AU ~ l/|Af/|). The position of the 
domain wall /* is determined by the condition that, at stationarity, the traveling time 
of particles from the shock region to site L must equal to the traveling time of vacancies 
through the high density domain to site 1, otherwise the domain wall would move with 
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a finite velocity. Thus, for ^, ^ <^ L, we may write 

Introducing the fraction r of the system which is occupied by the high density domain 
in the hmit L ^ oo, i.e. r = hm^^oo ^*(-^^)/-^, Eq. (ITTI) yields: 

.^(:.,AWA.)-'^(lH-,ni|^/.n*±i;|)-\ (IS, 

As can be seen, unlike for the one-lane partially asymmetric simple exclusion process, r 
differs from 1/2 if u ^ v, which is the consequence of the breaking of the particle- hole 
symmetry. For large L, the stationary current, which is proportional to the inverse of 
traveling time, tends exponentially to zero with increasing system size: 

J(L) ~ e-^^(-^)^ = f 4^ • (19) 

\p{q + u)J 

If /3 > (and a > 0), the high density domain shrinks to the left boundary of the 
system; in other words, the low density domain extends over the entire system and a 
boundary layer forms at the left boundary. The width of the boundary layer increases 
with decreasing j3 but remains finite for finite j3, which means that r = 0. The expression 
of the current in Eq. (fT9|) is valid with r = 0. 



5.2. Marginal bias 

Next, the marginal case is considered, when q/p = u/v or, equivalently, AU = 0. From 
Eq. f|T8l) we obtain r = 1. Therefore, if /? = and a > the high density domain 
extends over the whole system and a boundary layer forms at the right boundary. 
The characteristic scale ^ = 1/|A?7| diverges, therefore the density must tend slower 
than exponentially to the bulk value. The motion of a single vacancy is diffusive, 
so, far from the right boundary, where the density of vacancies is close to zero, they 
can be approximatively regarded as independent, diffusing particles. When a system of 
independent, diffusing particles is coupled to boundary reservoirs with different densities, 
the steady state density profile interpolates linearly between the densities of reservoirs. 
Therefore we conclude that, for sites far away from the right boundary, the profile is 
linear: 1 — (nf) i/L, 1 — (nf) ^ i/ i <^ L. It follows then that the current is 
inversely proportional to the system size: 

J{L) = a{l - {n^)) ~ L~\ (20) 

For /? > (and a > 0), the high density domain shrinks to the left boundary 
(r = 0). For a given /5, the width of the boundary layer is significantly larger than in 
the case of strong bias since the profile in the high density domain approaches the bulk 
value much slower. The current decreases exponentially as given in Eq. ( fT9l) with r = 0. 

Numerical simulations are in agreement with the above results obtained for strong 
and marginal bias. 
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5.3. Weak bias 

Next, we turn to the case of weak bias in the bulk. According to numerical simulations, 
the open system has three phases, see the phase diagram in Fig. El If the exit rate 




Figure 6. Phase diagram of the open system with weak (reverse) bias. The bulk rates 
are p = 2/3, q = I, u = 3/4 and v = 1/4. The low density, high density and maximum 
current phase are denoted by LD,HD and MC, respectively. The boundaries of the 
high density phase have been calculated by numerical simulations. Along the dashed 
line given in Eq. (|23p . the semi-infinite system has a product measure steady state. 
The circles denote the exactly known points of the phase boundaries. 



exceeds a certain a-dependent value, l3i{a), the effective injection rate of particles at 
the left boundary is not sufficiently large to force a non-vanishing current against the 
bias in the bulk. The low density domain covers the entire system (r = 0) and a 
boundary layer forms at the left boundary. The current decreases exponentially with L 
as given in Eq. (|T9l) with r = 0. This phase is termed as low density (LD) phase. If j3 
is smaller than l3i{a) the particle injection at the left boundary is able to force a finite 
density ^ > to the bulk. At both ends of the system, boundary layers form but, far 
from the boundaries, the profiles tend exponentially to constant values p and tt in lane 
A and B, respectively, which fulfill relation (fT2l) . Thus, in the bulk of a large system, 
the influence of boundary layers can be neglected. To be precise, the steady state of a 
finite open subsystem in the middle of the system, [L/2 — a,L/2 + a], is of factorized 
form in the limit L ^ oo: P{{nf,nf}) = IlflV/2-a P^'H^ - p)^""^vr"f (1 - 7r)i-"?. As 
a consequence, the current is given by Eq. (fT3l) in the limit L — oo and, in accordance 
with g > Qo, it is non-vanishing: Jqo = liniL^oo J{L) > 0. In this phase, which is 
termed as high density (HD) phase, the bulk densities p and tt, as well as the current 
are non-trivial functions of a and (3. 

At the phase boundary between the LD and the HD phase (/? = Pi{a)), a similar 
phenomenon can be observed as at the coexistence line of the one-lane ASEP. The HD 
and the LD phase coexist here in the sense that a delocalized anti-shock develops in the 
interior of the system, which separates a HD domain on its left hand side from a LD 
domain on its right hand side. At the left end of the system, exponentially decaying 
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boundary layers develop in the profiles. In the bulk of the HD domain, the densities 
are po ^^^d ttq given in Eq. f|T6|l . which follows from that J^o = 0. In the bulk of the LD 
domain, the densities are zero: p = vr = 0. The anti-shock has a finite width and its 
center performs symmetric diffusion in the whole system. A symmetric random walk on 
a one-dimensional lattice of size L with reflective boundaries has a steady state with a 
uniform probability 1/L at each site. For a sharp domain wall and no boundary layer 
at the left end this would lead to profiles which interpolate hnearly between po(7ro) and 
zero in lane A{B). For a finite system this holds only approximatively for large L and far 
from the boundaries: (nf) ~ {l—i/L)pQ, {nf) ^ {l — i/L)TiQ, 1 <^ i <^ L. Nevertheless, 
on the macroscopic scale, i.e. introducing the rescaled variable x = i/L and performing 
the limit L oo, the profiles are given as: 

(n^) = (1 - x)po, (nf ) = (1 - x)7ro (21) 

for < X < 1. The dependence of the current on the system size L can be easily 
obtained. Measuring the current at the right end of the system, it is clear that particles 
have a significant chance to leave the system at the right boundary and contribute to the 
mean current only when the domain wall is found in the vicinity of the right boundary. 
The probability of this event is proportional to L~^, therefore the current in a finite 
system decreases with L as: 

J{L) ~ L-\ (22) 

In the HD phase, the current Joo is determined by the boundary rates at the left 
boundary and, for a fixed a, it increases with decreasing /?. For given bulk rates, 
however, the current which can flow in the bulk is bounded and has a maximal value 
<^max5 see the fundamental diagram in Fig. HI Therefore, if a is not too small and /3 is 
decreased the current attains its maximum at some P = /32(«) > 0. If /?2(«) > and 
P < /?2(«), the current is saturated, i.e. Joo = Jmax, and the bulk densities pmax and 
TTmax in lauc A and B, respectively, are given by the equation Joo(Pmax) = -^max and by 
Eq. (|T2l) . These densities, as well as the current are independent of the boundary rates 
a and (3. The density profiles tend algebraically to the bulk value with the distance / 
measured from the boundaries as \{nf) — p| ~ \{nf) — tt ~ /^-^/^ for 1 <^ I <^ L. 

The value of the decay exponent 1/2 is related to that the fundamental diagram is 
quadratic in the vicinity of the maximum [10|. After the analogous phase of the ASEP, 
this phase will be called maximum current (MC) phase. 

In general, we could not determine the phase boundaries /3i(a) and /52(a) exactly 
apart from certain special cases. For a given set of bulk rates, we have performed 
numerical simulations for different boundary rates a and (3 and estimated the phase 
boundaries by measuring the current in the steady state. The results are shown in Fig. 
[6l The system size and simulation time was L = 300 and t = 2 x 10®, respectively, in 
case of the upper boundary Pi{a), whereas L = 1000 and t = 5 x 10® in case of the 
lower boundary /?2(a)- 

A special case, where the upper boundary can be determined exactly is p = q. If 
vacancies(particles) in lane B are regarded as particles (vacancies), one can see that none 
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of the two lanes is singled out if a = /5. The rates at the left boundary thus force a high 
density domain to the system in which {nf)HD = 1 — {nf)HD holds. Here, the subscript 
HD refers to that the average is restricted to configurations where site i lies in the high 
density domain, i.e. i <^ / where / is the location of the domain wall. Thus, in the bulk 
of the high density domain p = 1 — vr and the corresponding current is zero. Therefore 
the phase boundary curve Pi{a) coincides in this special case with the diagonal a = p. 

Furthermore, there is a line in the phase diagram, where the semi-infinite system 
has a product measure steady state. It is easy to see that if the densities of particle 
reservoirs p = a/p and vr = 1 — f3/q are compatible in the sense that they fulfill Eq. 
(fT2l) . moreover, the system is in the HD phase, i.e. po < a/p < pmax, then any finite 
open subsystem [l,a] (a < oo) has a factorized steady state in the limit L ^ oo with 
homogeneous densities p and vr in lane A and B, respectively. Expressing the densities 
of boundary reservoirs with the boundary rates and substituting them in Eq. ( |T2|) . we 
obtain 



This curve is plotted in the phase diagram in Fig. El On this line, there is no boundary 
layer at the left boundary. The upper phase boundary is reached at a = ppo, i.e. 
PiiPPo) = Pfactippo) , whereas the lower phase boundary is reached at a = ppmax, i-e. 
P2{.PPmax) = /5fact(pPmax)- This providcs oue exactly known point of each phase boundary 
curve, which are shown in Fig. [6l 

6. Inhomogeneous model 

In the rest of this work, we shall investigate the steady state of the model in the case 
where the system is heterogeneous, which means that the jump rates are site-dependent. 
If the bulk rates are random variables, one can easily imagine that there may be both 
regions where particles feel a local forward bias and regions with a local reverse bias. 
Such a situation has been studied in the disordered one-lane partially asymmetric simple 
exclusion process [lOl [35l [36l EH [38l [39]. Note that, in the totally asymmetric process, 
such a situation cannot be reahzed. The key to the understanding of such heterogeneous 
models at a phenomenological level is the steady state of a one-barrier system, i.e. when 
a homogeneous region with a reverse bias is embedded in a homogeneous medium with 
a forward bias. 

6.1. One-barrier system 

Let us consider a periodic system of size C, which is composed of a homogeneous 
subsystem (called barrier) with a reverse bias {AU > 0), which contains sites i = 
1,2, ... ,L and the rest of the system (the environment), where there is a forward bias 
{AU < 0). The particle density g = N/{2C) is fixed, furthermore, we assume that 
1 ^ L = bC, where 6 <^ 1, thus the barrier is only a small part of the system so that a 




(23) 



Dynamics at barriers in bidirectional two-lane exclusion processes 



14 



single particle has a positive average velocity (in a finite system). Similar to the reverse 
bias regime, one must distinguish between two types of forward bias. If AU < and 
AU > 0, we speak of strong forward bias, while if AU < but AU < we speak of 
weak forward bias. Depending on whether the bias is strong or weak at the barrier and 
in the environment we have four different cases. 

Let us assume that the bias is strong in both parts of the system. In this case, 
almost all particles are found in a high density cluster of length I ~ N/2 = gC, where 
p = TT = 1, see Fig. [7^. The right edge of the HD domain (the anti-shock) lies in 
the barrier region and covers a fraction r of this region as given in Eq. (ITSl) . provided 
the number of particles is sufficiently large: N/2 ^ I > 2rL. Focusing on the barrier 
region, the steady state is similar to that of an open system with a > and (3 = 0. If, 
however, N/2 < 2rL the HD cluster is located symmetrically on the two sides of site 1 
and covers only a fraction r = 1/{2L) = N/{4L) of the barrier region. The stationary 
current decreases exponentially with L as given in Eq. (fT9|) . 




d) n V 



Figure 7. Density profiles in different one-barrier systems. Arrows indicate tfie 
direction of the bias and big(small) arrow heads correspond to strong(weak) bias. 

Next, we consider the case when the bias is still strong in the reverse bias region 
but it is weak in the environment. Now, a HD cluster forms in the forward bias region, 
the right edge of which is located at site 1, see Fig. [Tb- In the HD domain, the 
densities po and ttq are given by Eq. (fT2|) with the rates of the forward bias region. As 
particles detaching from the anti-shock region must overcome the empty barrier region, 
the stationary current is exponentially vanishing, as given in Eq. (fT9l) with r = 0. 

Let us consider the case when the bias is strong in the environment but it is weak in 
the barrier region. Reversing the directions left and right, this arrangement is identical to 
the previous one, except that the barrier region is limited. Nevertheless, until < 2LgQ 
or, equivalently, if the length / = N/{2qq) of a HD cluster with density qq is less than 
L, a HD domain will form in the reverse bias region as we have known in the previous 
case. This is illustrated in Fig. [T^. At the left edge of the HD domain (at site 1), a 
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boundary layer forms in the forward bias region, which effectively closes the reverse bias 
region, and the situation is similar to the open system with a = j3 = 0. Since, the right 
edge of the HD cluster does not reach the right edge of the barrier region, the current is 
exponentially vanishing as given in Eq. (fT9l) with r = l/L. The steady state is, however, 
different if > 2Lgo. In this case, the HD domain covers the entire reverse bias region 
and, here, the density exceeds Qq. Besides, the density of particles is no n- vanishing also 
in the environment and, accordingly, the current is non- vanishing (Jqo > 0), as well. 

The most complex situation is when the bias is weak both in the barrier region and 
in the environment. Restricting our attention to the reverse bias region, this can be 
regarded as an open system of size L, with some effective boundary rates Ocs and (3ea, 
which depend on the bulk rates in a non-obvious way owing to the presence of boundary 
layers around site 1. Since both effective boundary rates are finite, the determination 
of the steady state in the general case is, in addition to the above difficulties, as hard 
as the determination of the phase diagram for an open system with weak reverse bias. 
We shall therefore restrict ourselves to the description of the possible steady states and 
omit to give criteria in terms of the jump rates except of some special cases when the 
symmetries of the model can be exploited. 

For almost all points of the parameter space spanned by pf,qf,Uf,Vf,Pr,qr,Ur,Vr, 
where the index /(r) refers to the forward(reverse) bias region, the steady state is 
similar to those described in the previous two cases. That means, a high density cluster 
forms which is located either entirely in the forward bias region (Fig. [Tb) or entirely in 
the barrier region (Fig. [7b) and the densities are given by the densities po and ttq of the 
medium which correspond to zero current, see Eq. (fT6|) . In the latter case, the current 
can be non-vanishing (Joo > 0) if the density of particles is sufficiently high. 

In addition to the above steady states, there is a zero measure set in the parameter 
space for which the high density cluster becomes delocalized, analogous to the behavior 
at the phase boundary between the LD and the HD phase in the open system (see Fig. 
[TJi). In this case, the right edge of the high density cluster lies in the reverse bias region 
and the whole cluster performs symmetric diffusive motion so that the right edge of 
the cluster moves in the range [1, min{L, Z^}], where Ir = N/{2qq) is the length of the 
cluster corresponding to the zero-current density in the reverse bias region. In the 
part of the cluster which lies in the forward (reverse) bias region, the local densities 
assume the zero-current densities of the underlying medium, i.e. pi and ttq (pq and tTq). 
If Ir/L = f?/(6f?o) — then the current vanishes with the size of the barrier region as 
J{L) ~ L^^, otherwise it decreases exponentially. The symmetric diffusive motion of 
the cluster results in a linearly decaying average density profile in the barrier region. 

The general condition of such a delocalized steady state is not known but we can 
provide two sufficient conditions, which are related to the symmetries of the model. 
First, if pf = qr, qf = Pr, Uf = Vr and Vf = Ur then the two sides of site 1 (i.e. the 
forward bias region and the reverse bias region) are mapped to each other by rotating 
the system by 180 degrees about the axis perpendicular to the plane in which it is drawn 
in Fig. [H That means that the two regions are equivalent and the currents in the two 
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lanes cancel, independent of the actual location of the HD cluster. Fluctuations of the 
net current lead to displacements of the latter and, as a consequence, it is delocalized 
in the steady state. Second, if Pf = qj = Pr = Qr, the high density cluster is delocalized 
again, independent of the lane change rates. This latter case is related to the symmetry 
already mentioned in Sec. 15.31 When the particles(vacancies) in lane A are regarded 
as vacancies(particles) then the particles in both lanes move in the same direction 
and lane change of the original model transforms to pair annihilation and creation. 
If the above condition is satisfied the lanes in the transformed model are equivalent. 
Therefore, the particle currents in the two lanes of the original model compensate each 
other independent of the location of the particle cluster, which leads to delocalization. 
Although, boundary layers develop around site 1 in the profiles in lane A and B, for the 
total density {nf)HD + {nf)HD = 1 holds in the HD cluster due to the above symmetry. 

We have performed numerical simulations in the second case and measured the 
stationary density profiles (see Fig. [8]), as well as the current for different system 
sizes. The results on the current, which are shown in Fig. [9l are in agreement with 
the law J{L) ~ L^^ for large system sizes. The intercept L* of the linear asymptote 
J^^[L) ~ s{L — L*) with the x axis can be interpreted as the effective width of the 
anti-shock region. As can be seen, for smaller Vr/ur, the anti-shock is sharper and the 
finite size corrections to the asymptotic law are weaker. The slope s of the asymptote 
depends on Vr/ur, as well, through the sharpness of the anti-shock. 
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Figure 8. Numerically calculated density profiles in a periodic one-barrier system. 
The size of the system is £ = 800, the number of particles is = 100 and 
the reverse bias region is located in the interval [200,300]. Transition rates are 
Pf — Qf = Pr = ^ 0.5 and Vr = 1 — Ur = Uf = 1 — Vf. The simulation time 
was 10^. Only the part of the system is shown where the density is significantly differs 
from zero. In the rest of the system, which is never reached by the high density cluster, 
the density is vanishing: p « J/iPfi'^^f ^ C~^. 

6.2. Disordered system 

Next, we shall study the model in the presence of disorder, where the transition rates 
are site-dependent quenched random variables. The general case is rather complex 
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Figure 9. The inverse of the measured stationary current as a function of the size L 
of the barrier in a periodic one-barrier system for different lane change rates. The size 
of the system was C — 16L and the number of particles was N = 4L. Transition rates 
are Pf = q/ = Pr = Ir- = 0.5 and Vr = 1 — Ur = Uf = 1 — Vf. The simulation time was 
109. 



as, in such systems, both weak and strong bias regions may form. We shall therefore 
restrict the quantitative analysis to two special cases where either exclusively strong or 
exclusively weak bias regions are found in the system. 

6.2.1. Longitudinal disorder First, let us consider the case where the lane change rates 
are homogeneous, i.e. Ui = Vi = u {i = 1,2 L) while the other rates pi and qi are 
independent quenched random variables. The model with such a longitudinal disorder 
may describe traffic of molecular motors on heterogeneous tracks containing defects 
where the advance of motors is hindered. 

In this case, the potential difference AUij between site i and j(> i) felt by a particle 
and the corresponding potential difference AU ij of a vacancy are related as: 



The second term on the r.h.s. is an 0(1) random variable (with zero mean), which is 
negligible for large j — i compared to the first term which increases linearly with j — i 
in the case of a global bias and increases as ~ y/j — i in the absence of a global bias. 
Thus the two potentials are asymptotically perfectly correlated: 



As a consequence, there are only strong (forward or reverse) bias regions but no weak 
bias regions in the case of longitudinal disorder. This is the case also in the disordered 
one-lane partially asymmetric simple exclusion process, therefore the model behaves 
similar to that one and we can make us of the known results on that model. 

In the following, we assume that the potential landscape is tilted, i.e. [Af/jjav < 0, 
where [-Jav denotes average over the distribution of jump rates. A given realization of 
the random potential U^, i = 1,2, ... ,L can be regarded as a random walk of L steps 




(24) 




(25) 
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with step lengths AUi, in the presence of a bias in the negative direction. The effective 
reverse bias regions in the potential landscape can be identified with the ascending parts 
of Brownian excursions in the positive direction [371139]. A Brownian excursion of length 
n is a random walk Ui, i = 0,1, ... ,n, which satisfies that = Uq and f/, > Uq for 
< i < n. The amplitude of the excursion is defined as f/max — Uq = maxo<i<n Ui — Uq. 
The effective reverse bias region is identified with the part of the excursion from i = 
up to the first maximum: 

/= min{j|f/, = [/^ax}. (26) 

0<j<n 

In the steady state of a finite system, the current is determined by the reverse 
bias region with the largest amplitude AUmax, behind which the particles accumulate 
and form a high density cluster, where the density is close to 1. Due to the particle- 
hole symmetry, which is reflected by Eq. fl25l) . half- filling is realized in this maximal 
reverse bias region, i.e. the anti-shock is located where the potential (measured from the 
bottom of the barrier) is A[/max/2. The current is therefore given by J = I/^tmax, 
where tmax = e^^'^^'^ is the waiting time of a single particle in the maximal reverse bias 
region. The distribution of the amplitudes of Brownian excursions decays asymptotically 
exponentially (see the Appendix) or, equivalently, the probability density of one-particle 
waiting times at reverse bias regions has an algebraic tail p(r) ~ At~^~'^. Here, the 
diffusion exponent fi is the positive root of the equation [(g^/p)^]av = 1 in case of a 
one-lane system with independent, random rates pi and (ji (see e.g. Ref. [H]). In terms 
of the jump rates of the original two-lane model, this can be written as: 
' q{p + u) \ 
^p{q + v)) 

According to extreme value statistical considerations [12], the current, which is related 
to the largest one-particle waiting time, follows the Frechet distribution 

p{J) = 2/ij2A«-ie-^'"*, (28) 

in terms of the scaling variable J = cJ{L)L^ , where the constant c is related to the pre- 
factor A of the parent distribution p{t). The typical current vanishes with the system 
size L as: 

J{L) ^ L-^/^^f^K (29) 

We have performed numerical simulations for the model with rates u = 0.5, Pi = 1 — qi 
for all i and rate pi was drawn from a bimodal distribution 

p{p) = c6{p - A(l + A)-^) + (1 - c)6{p - (1 + A)-^), (30) 

with A = 0.25 and c = 0.2. For this bimodal randomness, the exponent /i is given by 

X+{1 + X)u 



1. (27) 



ln(c 



In 



(31) 



A + A(1 + A)m_ 

which yields l/(2yu) = 0.2767288 . . . for A = 0.25 and c = 0.2. We have measured the 
stationary current in 10^ Monte Carlo steps in 10^ independent samples for system sizes 
L = 2*, 2 = 8, 9, 10, 11, 12. The results on the distribution of the current shown in Fig. 
[TUl are in agreement with the above predictions. 
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Figure 10. Left: Distribution of the logarithm of the current obtained by numerical 
simulations for different system sizes in case of the bimodal longitudinal disorder 
given in Eq. ((30)) with A — 0.25 and c = 0.2. Right: Scaling plot with the 
exponent l/(2/j,) = 0.2767288 • • • predicted by the theory. The solid line is the Frechet 
distribution given in Eq. (|28|) . 



6.2.2. Transversal disorder Let us now consider another type of randomness, where 
the longitudinal rates are homogeneous, i.e. pi = qi = p for all i, whereas the transversal 
rates Ui and Vi are random variables. This model will be termed transversal disorder and 
may describe situations where the lane change of molecular motors is biased randomly 
in space due to the presence of forces transversal to the filament. 

In this model, the one-particle potential and one-vacancy potential are again 
perfectly correlated: 

AUi = AUi. (32) 

Consequently, there are exclusively weak (forward or reverse) bias regions in the 
system, while strong bias regions are lacking. Let us assume in the following that 
the potential landscape is tilted, i.e. [Af/j]av < 0. Similar to the longitudinal disorder, 
the environment in which the particles move can be regarded as a set of isolated reverse 
bias regions which are embedded in a medium with forward bias. We shall apply a 
phenomenological independent-barrier model to this system, some elements of which 
have already been mentioned in the previous section. This theory which gives exact 
results for a single particle in a random environment and results which are in agreement 
with numerical results for the disordered partially asymmetric simple exclusion process 
|37] . as well as for the present model with longitudinal disorder. The system will be 
represented by a set of L' = 0{L) consecutive barriers, which exchange particles with 
each others. There are A^, particles at the ith barrier and the barriers are characterized 
by own current functions /j(iVj) which give the current in the one-barrier system with 
Ni particles, where the other barriers are deleted. The current /j(A^j) is monotonically 
increasing with Ni up to a saturation value /™^^, i.e. Ii{Ni) = If^^^ if Ni is greater than 
some N^^^ which is characteristic of the barrier. In what follows, we assume that the 
total number of particles is sufficiently large meaning that X^i Nf^^^ < N. Then, in the 
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steady state, the current is determined the barrier with the minimal carrying capacity 
in the system (indexed by m): J{L') = minjjJ™'*^} = I^^^. The mean number of 
particles Ni at other barriers is determined by Ii{Ni) = J{L') and Ni < N^^^ for i ^ m, 
whereas the number of particles at the saturated barrier is Nm = N — J2i^m > ^m^^- 
These particles form a high density cluster behind the saturated barrier, see the density 
profiles in Fig. [TTJ As we have seen in the previous section, the transversal disorder 




100 200 
i 

Figure 11. Numerically calculated density profiles in the steady state in the presence 
of transversal disorder. The system size is L = 200 and the number of particles is 
N = 100. Bimodal randomness has been used given in Eq. (|33p with c — 1/3, A = 1/8 
and p = 0.5. The dotted line indicates the potential landscape Ui. The longest reverse 
bias region is in the domain [139, 170]. Behind this region, a high density cluster of 
approximate length 100 can be observed, where (nf) + {nf) w 1. 

corresponds to the special case, where the high density cluster is delocalized at a barrier. 
For homogeneous barrier regions, which is a good approximation for bimodal randomness 

of the type = 1 — t>j, 

p{u) = c6{u - A(l + A)-i) + (1 - c)6{u - (1 + A)-i) (33) 

in the limit c <^ 1, the saturation value of the current (when Ni > Nf^^^) depends 
asymptotically on the extension of the barrier region I as [/™'^^(Z)]~^ ~ s{l — /*) where 
s and /* are common constants for the barriers. Apart from this limit, the barriers 
are inhomogeneous, i.e. they may contain smaller domains with forward bias, and as a 
first guess, they are identified with the ascending parts of Brownian excursions of the 
potential Ui, just as for longitudinal disorder. Nevertheless, for inhomogeneous barriers, 
the HD cluster is still delocalized, therefore the saturation value of the current scales 
with the extension li of the barrier defined in Eq. ( l26l) as I^^^(li) ^ Thus, the 
current in a finite system is expected to scale as 

J - Cx, (34) 
where /max = niaxi<j<i/{Zj} is the length of the longest reverse bias region present in 
the system. 
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We have performed numerical simulations for the model with transversal disorder 
and measured the current in the steady state. We have determined in each sample the 
size /max of the longest effective reverse bias region. Then the measurement was carried 
out for 10^ independent samples and the average of the logarithm of the current has 
been calculated over samples where /max was the same. The dependence of the calculated 
typical current [J(/max)]typ = exp([ln J(/max)]av) on /max is shown in Fig. [121 As can be 




'max 



Figure 12. Dependence of the typical current on the length of the longest reverse 
bias region in the model with transversal disorder. The size of the system was L = 256 
and bimodal randomness in Eq. ([55)) has been used with c = 1/3 and A = 1/8. The 
number of particles was N — 128. 

seen, the numerical results are in agreement with the law given in Eq. (IMl) . Independent 
of the precise definition of effective reverse bias regions, one may generally conclude that 
their probability of occurrence must be exponentially small in their size / (apart from 
possible corrections like in the case of continuous-time Brownian excursions [Hj) since 
the majority of the lane change rates in such a region is of that type which results in 
reverse bias (in the particular case in Eq. ([33|1 that with probability c). Therefore, the 
length of the longest effective reverse bias region scales with the size of the system L as 
/max ~ In L. Finally, we obtain that the current scales with L in the case of transversal 
disorder as 

J{L)^{\nL)-\ (35) 

We have computed the probability of amplitude AU of Brownian excursions and the 
probability of the extension / of the effective reverse bias region which is given by the 
time at which the maximum is first assumed in Brownian excursions; for the details of 
the calculation, see the Appendix. The results are plotted in Fig. [131 The probability 
of amplitudes, which is the relevant variable in case of longitudinal disorder, tends 
rapidly to the asymptotical exponential from. Contrary to this, the probability of the 
extension / of effective reverse bias regions, which is the important variable in case of 
transversal disorder, has strong corrections to the exponential for moderate /. Therefore 
the asymptotic decay in Eq. ([35ll is expected to be observed for rather large system 
sizes which are beyond the realm of numerical simulations. 
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Figure 13. Probability of the amplitude AU (in units In p^i/[x^^) ) and the 
probability of the time at which the maximum is first assumed in Brownian excursions 
corresponding to bimodal randomness with c = 1/3. Latter quantity is identified with 
the length / of effective reverse bias regions. The straight line is the asymptotic form 
given in Eq. IA.3I 



7. Summary 

We have studied in this work a two-lane model that describes bidirectional transport 
of particles interacting by hard core exclusion. Although, the one-particle dynamics in 
such a two-channel environment are qualitatively similar to those in a one-lane partially 
asymmetric system which allows backward steps of the particle, the many particle 
steady-state behavior can be much different in the two models. We have investigated 
situations when the track contains regions where the preferred direction of motion of 
particles is opposite to the global bias. In this case, we have found two types of steady 
states and have given the conditions of them in terms of the transition rates. In the 
case of strong reverse bias, which can be induced by making the rates of transitions 
parallel to the track asymmetric, the steady state is similar to that of the one-lane 
partially asymmetric exclusion process: a compact, localized cluster of particles forms 
and the current vanishes exponentially with the extension of the reverse bias region. In 
case of weak bias, realized by rendering the lane change rates asymmetric, which has 
no counterpart in the one-lane asymmetric model, a qualitatively different steady state 
is observed. The cluster of particles becomes now delocalized and, as a consequence, 
the current vanishes inversely proportionally with the size of the barrier region. We 
have thus found that, compared to analogous one-lane systems, there may be physically 
different steady states, in which the flow against the local external drive is facilitated 
by the cooperative behavior of particles. It is an interesting question whether this type 
of delocalization can be observed in experiments or in living cells. 
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Let us consider tlie bimodal randomness given in Eq. (!30l) and introduce Xi = 
p+i/{i+A) • "^^^^ rescaled potential landscape {xi} is a random walk with steps 
of unit length with probabilities c and 1 — c in the positive and negative directions, 
respectively. Let us assume that the walk starts at a;o = 1 and absorbing walls are 
located at x = and atx = X + l>l. A central quantity in the computation of the 
statistics of the maximum is the probability that the walker is found after t steps at site 
x: Px{x,t] c). The probability that the walk ends at the absorbing site X + 1 without 
ever having crossed the starting point is called persistence probability and it is given as 
Ppers(-^;c) = limt^oo Px{^ + ^,t',c). This quantity can be exactly calculated [13] and 
reads in the present case: 

Ppe.s(X;c) = (1/c- 2)[(l/c- 1)^+^ - (Al) 

The probability that the amplitude x^ax of the walk before it crosses the wall at a; = 
can be written in terms of the persistence probability as 

Prob(xmax = x) = c[ppcTs{x - 1; c) - Ppers{x; c)] (A. 2) 

for X > 1. For x = and x = 1, we have Prob(xmax = 0) = 1 — c and 
Prob(xmax = 1) = c(l — c). For large x, this probability decays exponentially: 

Prob(x^ax = x) ^ (1 - 2c)V(l - c)(l/c - 1)-^ (A.3) 

Next, we turn to the computation of the distribution of the time tmax at which the 
walk first reaches its maximum before it crosses the wall at x = 0. It is easy to see 
that Prob(tmax = 0) = 1 — c and Prob(tinax = 1) = c(l — c). For tmax > 1, the joint 
probability that the amplitude is Xmax and it is assumed for the first time at time tmax 
is given as c^P^n^ax-ila^max - I, tmax - 2; c)ppers(a;max; 1 - c). From this, the marginal 
distribution of tmax is obtained as 

t 

Prob(tmax = = J2 ^xn,ax-l(a;max-l,t-2;c)ppers(a;max; 1-C).(A.4) 

The quantities Px{x,t;c) which appear in the above expression can be computed 
iteratively from those at time t — 1 using the initial condition Px{l,0;c) = 1 and 
Px{x,0]c) = for X ^ 1. 
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